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Density-functional study of the pure and palladium doped small copper and silver 

clusters 
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The size-dependent electronic, structural, magnetic and vibrational properties of small pure cop¬ 
per and silver clusters and their alloys with one and two palladium atoms are studied by using 
full-potential all-electron density functional computations. The stable isomers of these clusters are 
identified and their theoretical magic numbers are determined via the analysis of the second differ¬ 
ence of their minimized energy. We discuss that the doped Pd atoms generally prefer to sit in the 
high coordination sites of the pure clusters. It is argued that Pd doping influences the structural 
properties and the two dimensional to three dimensional structural cross over in the small Cu and 
Ag clusters. The many body based GW correction is applied for more accurate determination of 
the electron affinity and ionization potential of these systems. Magnetic and vibrational properties 
of the pure and doped clusters are presented and discussed. 


I. INTRODUCTION 

Noble metal nanoclusters are currently attracting a 
great deal of attention within both experimental and 
theoretical communities, due to their interesting struc¬ 
tural and electronic properties and promising technologi¬ 
cal applications in nano electronics, nano optics, biologi¬ 
cal sensing, catalysis, and biomedicine [H-Q. While tran¬ 
sition metal ions may induce health and environmental 
problems SH , noble metal nanoclusters exhibit lower 
toxicity and hence more biocompatibility. Noble metal 
nanoclusters are alloyed with palladium atoms to im¬ 
prove their catalytic activity [1,1^. In recent years, ex¬ 
tensive computational studies have been devoted to prop¬ 
erties of pure noble metal clusters [IMl, although less 
information is available about structural and electronic 
properties of the corresponding bimetallics [l^-[T^ . 

Among noble metal based bimetallic nanoclusters, cop¬ 
per and silver based nano-alloys have been studied con¬ 
siderably less than gold based bimetallic nanoclusters. 
Wang et al. El and Romanowski et al.Ell used den¬ 
sity functional (DFT) computations to study AgPd and 
CuPd dimers and their interaction with H 2 . Other cal¬ 
culations have been performed on Cu-Pd trimers and 
their interaction with molecular and atomic oxygen E3- 
The stable geometry of small Ag„Pdm clusters up to five 
atoms were found to transform from two-dimensional to 
three-dimensional as the Pd content increases [l^. Zhao 
et al. performed first-principles calculations to study the 
effects of Pd doping on structural and electronic prop¬ 
erties of Ag„Pd (n < 5) clusters and their mono hy¬ 
drides [2l| . Efremenko et al.have studied the geomet¬ 
ric structures and electronic properties of the Pd„Cum 
(n -I- TO < 6) clusters using DFT. The obtained results 
show that stability of nanoclusters of the same shape and 
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composition increases linearly with increasing number of 
PdCu bonds. Up to our knowledge, there is no report 
on magic number, vibrational spectra, and many body 
corrected electronic structure of small Pd doped Ag and 
Cu clusters. 

In this paper, quantum mechanical calculations are 
employed to study pure M„ (M = Ag, Cu) and doped 
M„_mPdm {n < 9,m = 1,2) clusters and investigate 
the behavior of their structural, electronic, and vibra¬ 
tional properties as a function of size. The rest of the 
paper is organized as follows. In section II, we give a 
brief introduction of the computational method used in 
this work. Then in section III, atomic dimers are inves¬ 
tigated with several exchange-correlation functional to 
select the proper one for our calculations. Next section 
involves the results of our search for lower energy struc¬ 
tures of the selected nano-clusters. In section IV, the 
structural, electronic, magnetic and vibrational proper¬ 
ties of the most stable clusters are presented. Finally, we 
will summarized our findings in the last section. 


II. METHOD 

All presented calculations are performed in the frame¬ 
work of the spin-polarized Kohn-Sham density functional 
theory by using the all electron full-potential code FHI- 
AIMS Q. This package employs basis sets consisting of 
atom-centered numerical orbitals of the form: 

«(b = (1) 

r 

As the name implies, the radial free atom like orbitals 
Ui{r) are numerically tabulated and therefore exhibit 
yery high flexibility. The strict localized nature of these 
basis functions lead to much slower scaling of the compu¬ 
tational time versus the system size. The calculations re¬ 
ported here are done with the ”tier2-|-spd”, ”tierl-|-spd”, 
and ”tier2” basis sets which contain 116, 152, and 45 ba¬ 
sis functions for Cu, Ag, and Pd atoms, respectively. All 
the calculations were carried out in the scalar relativistic 
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limit, while the spin-orbit correction was neglected. A 
recent computational study on the atomic structure of 
transition metal clusters concluded that this relativistic 
term has negligible effect on the relative energy differ¬ 
ences of the 4d systems [24 1. 

Geometry optimization is performed by the stan¬ 
dard Broyden-Fletcher-Goldfarb-Shanno (BFGSj algo¬ 
rithm [ 2 ^ with a force accuracy of about 10“^eF/A. Har¬ 
monic frequencies are calculated using the finite displace¬ 
ment of all atomic positions by 10“^A. The lowest energy 
structures are confirmed to be the true minima by calcu¬ 
lating their vibrational frequencies. The FHI-aims pack¬ 
age enables us to describe electronic single-quasiparticle 
excitations in molecules by using many-body correction 
GW self-energy [ 2 ^. 


III. DIMERS 

First we focus on the atomic dimers to select a reli¬ 
able exchange-correlation (XC) functional for our com¬ 
putations and moreover to gain some insights about in¬ 
teratomic interaction in our systems. The equilibrium 
bond length, binding energy, and vibrational frequency 
of the Gu 2 and Ag 2 dimers were calculated W using dif¬ 
ferent XG approximations, including PBE [^, revPBE 
[ 2 ^ and BLYP [ 2 ^ generalized gradient (GGA) func¬ 
tionals and PW local density functional. The ob¬ 
tained results are compared with available experimental 
data in table ID It is seen that PBE, revPBE, and BLYP 
give significantly more accurate binding energy and vi¬ 
brational frequency while accuracy of equilibrium bond 
length is qualitatively the same within LDA and GGA 
based functionals. In the case of Gu 2 , the overall accu¬ 
racy of revPBE seems to be slightly better than PBE, 
while the parameters of Ag 2 within PBE is overall closer 
to experiment. Although, BLYP gives better binding en¬ 
ergy for Ag dimer, for larger silver clusters, other theo¬ 
retical calculations [3l| argue that the BLYP functional 
fails to predict the correct stable structure of the system. 
Therefore, we adapt the PBE and revPBE functionals for 
calculation of Ag and Gu based clusters, respectively. 

For better understanding of the interatomic bonds in 
our target systems, we calculated the Pd 2 , GuPd, and 
AgPd dimers and presented their results in table ID It is 
observed that the bond length is increasing in the Gu 2 - 
GuPd - Pd 2 series while the absolute binding energy and 
vibrational frequency is decreasing. Hence, Gu-Cu bond 
has a larger strength and stiffness compared with the 
Gu-Pd and Pd-Pd bonds. Bond length and normal mode 
frequencies of AgPd and Ag 2 dimers are almost equal 
while these are larger than those of Pd 2 dimer. The larger 
bond length of Ag 2 may be due to the large atomic radius 
of Ag relative to Pd and Gu atoms. Binding energy of 
Pd 2 is smaller than that of AgPd which in turn is smaller 
than Ag 2 dimer. 

The Mulliken population analysis of valence orbitals 
shows that Gu and Ag have a similar electronic configura- 


TABLE I: Calculated equilibrium bond length Te (A), binding 
energy Eb (eV/atom) and harmonic vibrational frequency uj 
(cm“^) of the Cu 2 , Ag 2 , Pd 2 , CuPd, and AgPd dimers. The 
Cu 2 and Ag 2 dimers are calculated in four different exchange- 
correlation (XC) functionals. The corresponding experimen¬ 
tal values (Expt.) are also given for comparison. 



XC 

Te 

Eb 

UJ 

CU2 

PBE 

2.21 

-1.13 

269 


BLYP 

2.23 

-1.11 

261 


revPBE 

2.24 

-1.02 

259 


PW 

2.15 

-1.39 

299 


Expt. 

2.21 [3^ 

-1.04 [33] 

265 [33] 

Ag2 

PBE 

2.57 

-0.91 

180 


BLYP 

2.60 

-0.83 

172 


revPBE 

2.60 

-0.79 

171 


PW 

2.48 

-1.45 

270 


Expt. 

2.53 [3^ 

-0.83 [3^ 

192 [^ 

CuPd 

revPBE 

2.32 

-0.79 

232 

AgPd 

PBE 

2.56 

-0.70 

179 

Pd2 

PBE 

2.49 

-0.49 

195 


Expt. 

2.48 [3^ 

-0.51 [37] 

210 [38] 


tion of gO.94^0.05^10 .01 dimers. Com¬ 

paring this configuration with the free Gu and Ag atoms 
electronic configuration evidences a small spd 

hybridization in these dimers. On the other hand, com¬ 
paring the electronic configuration of Pd in the Pd 2 dimer 
(g0.58p0.05^9.37) free Pd atom indicates 

a significant spd hybridization in the palladium dimer. 
The reason is that the valence shell of the free Pd atom 
is composed of fully occupied and fully unoccupied or¬ 
bitals and hence interatomic bonding between palladium 
atoms requires significant promotion of electrons from 
the occupied d states to the unoccupied s shell. While, 
the valence shell of the free Gu and Ag atoms has a half 
filled s orbital which is used for interatomic bonding. It 
is also found that the charge transfer in the Pd-Gu and 
Pd-Ag bonds happens from Gu and Ag to the more elec¬ 
tronegative Pd atom. 


IV. STABLE ISOMERS 

In order to investigate properties of atomic clusters, 
the first essential step is identification of the lowest en¬ 
ergy structure of the clusters. Therefore, we performed 
a careful search for the stable structure of the pure and 
doped clusters. In the case of the pure clusters, all prob¬ 
able atomic configurations of the clusters were included 
in our search. After accurate atomic relaxation of the 
relevant configurations and comparing their minimized 
total energies, the most stable isomers of the pure Gu„ 
and Agji clusters (n < 9) were identified and sketched 
in Fig. [TJ These findings are in agreement with a re- 
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FIG. 1: Obtained lowest energy isomers of the pure Ag„ 
and Cu„ and doped Cu„_iPd, Ag„_iPd, Cu„_iPd 2 , and 
Ag„_iPd 2 clusters for 3 < n < 9. The Pd atoms are shown 
by light balls while the Cu/Ag atoms are indicated by black 
balls. 


cent study on silver clusters [s^. Because of the pres¬ 
ence of two different atoms, the process of finding the 
global minimum energy structure of the doped clusters 
is more complicated than the pure clusters. Therefore, 
for finding the most stable structures of the palladium 
doped clusters, we tried to limit our search to the more 
stable configurations proposed by previous studies on the 
Cu„Pdm (n+m < and Ag„Pdm (n-|-m < 5l [^l40j| 
clusters. The resulting most stable isomer of the doped 
clusters are presented in Fig. [T] It is seen that the most 
stable M„_ 2 Pd 2 clusters are made of the maximum num¬ 
ber of pyramids and moreover the two Pd atoms in these 
systems tend to bond together. These observations were 
used for obtaining the most stable structure of the M 6 Pd 2 
and M 7 Pd 2 clusters. 

The results show that up to the size of 6 (heptamer), 
pure clusters prefer 2 dimensional (2D) planar structures 
while larger clusters stabilize in 3D geometries. In or¬ 
der to understand the origin of this 2D-3D cross over, 
we analyzed the average bond length (d) and the average 
coordination number of atoms {tic) in the most stable 
2D and 3D isomers of the Cu„ and Ag„ clusters (Fig.[2). 
It is seen that the 2D isomers have lower average coor¬ 
dination and shorter average bond length, because the 
valence charge density of these systems is distributed 
among lower number of bonds, compared with the 3D 
isomers. As a result of that, individual bonds in the 
planner isomers are stronger and shorter, while 3D iso¬ 
mers have more number of bonds. The 2D-3D structural 
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FIG. 2: Average bond length d and coordination number of 
atoms Tie in the most stable 2D and 3D structures of Cu„ 
(left) and Ag„ (right) clusters. 


transformation reflects the competition between the in¬ 
dividual bond strength and the total number of bonds. 
It is observed that (Fig. [2]) in the smaller sizes (n < 6), 
the more pronounced difference in the individual bond 
strength leads to the stability of the 2D isomers, while 
for the larger clusters the difference in the average co¬ 
ordination wins the competition and stabilizes the 3D 
isomers. 

The results (Fig. [ij indicate that the most stable 2D 
isomers of Cu„ and Ag„ clusters are made of equilateral 
triangles. It is also seen that the 2D-3D cross over in the 
Ag„_iPd clusters occurs at n = 6 while the Cu„_iPd 
clusters, follow the same structural transformation as the 
pure Cu cluster. As a general trend, we observe that 
the Pd atom prefers the high coordinated sites of the 
Cu„_iPd and Ag„_iPd clusters. The obtained results 
indicate high tendency of the M„_ 2 Pd 2 clusters toward 
3D configurations. Considering various 2D and 3D struc¬ 
tures for the M 2 Pd 2 , it was found that the 3D structures 
have lower energies, in agreement with the previous works 
[13, 113 . Based on these results, the theoretical search for 
the stable isomers of the M„_ 2 Pd 2 clusters were limited 
to all 3D isomers of the pure Ag and Cu clusters. It is 
worth while to mention that other authors have not in¬ 
vestigated Cu„_ 2 Pd 2 and Ag„_ 2 Pd 2 clusters larger than 
five and six atoms, respectively. 


V. STRUCTURAL PROPERTIES 

In order to address stability and structural behaviour 
of the clusters, the binding energy per atom, average co¬ 
ordination number of atoms, and the average bond length 
of the most stable structures of all clusters under study 
as a function of the cluster size are presented in Fig. 
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|3l The harmonic vibrational frequencies of the clusters 
were also calculated to address dynamical stability and 
bond stiffness of the systems. The maximum vibrational 
frequency of the clusters are shown in Fig. [31 The cal¬ 
culated binding energy of the copper clusters are in well 
agreement with the measured data. It can be seen that 
by increasing the cluster size, the binding energy of the 
pure clusters increases gradually toward the binding en¬ 
ergy of bulk copper (3.50 eV)[4|| and silver (2.95 eV)[4l|, 
although these clusters are clearly far from the bulk limit. 
Notably, the eight atom pure clusters occurs in a local 
maximum of the absolute binding energy. It is clearly 
seen that the 2D to 3D structural transformation in the 
pure and mono-doped clusters is reasonably accompanied 
by a sharp increases in the average coordination number 
of atoms. 

Observation of no imaginary frequency indicates the 
dynamical stability of the lowest energy isomers shown 
in Fig. |T] It is generally seen that the Ag based clusters 
have lower vibrational frequency than Cu based clusters, 
which is clearly related to the higher atomic mass of sil¬ 
ver. Considering the trend of maximum frequencies along 
with the geometry of the most stable isomers of the clus¬ 
ters (Fig. [T]) clarifies that the structural cross over from 
planner to 3D geometries in the pure Ag and Cu clus¬ 
ters and doped Cu„_iPd cluster is accompanied with a 
significant mode softening in the vibrational spectra. 

Comparing the properties of the pure and doped clus¬ 
ters show that Pd doping in the copper clusters increases 
the average bond length of the system, because of the 
larger atomic radius of Pd compared with Cu. This in¬ 
crease is expected to be accompanied by a bond softening 
effect, which is generally visible in the calculated har¬ 
monic vibrational frequency of the doped Cu clusters. In 
contrast to this bond elongating and softening effects, we 
observe that Pd doping has a week effect on the binding 
energy of Cu clusters. It seems that the doped Pd atom in 
the Cu cluster, inserts effective amount of electrons in the 
neighboring bonds and hence compensate the observed 
bond softening effects. It is reasonable, because Pd has 
much more electrons than Cu. On the other hand, the 
average bond length of the doped Ag clusters is slightly 
smaller than those of the pure clusters. The reason is 
that atomic size of silver is larger than that of palladium 
and hence Ag-Pd and Pd-Pd bonds are shorter than Ag- 
Ag one (table HI). As a result of that, a bond hardening 
effect is visible in the calculated vibrational spectra of 
the doped Ag clusters, compared with the pure ones. As 
a result of these bond shortening and hardening effects, 
Pd doping enhances the absolute binding energy of Ag 
clusters (Fig. [3]). 

In cluster science, the second-order energy difference 
{A 2 E) is a conventional quantity to address the relative 
stability of clusters as a function of size and reproduce the 
experimental magic numbers observed in mass spectra 
measurements. Large values of A 2 E{n) indicate that the 
n-atoms cluster is more stable than n — 1 and n-l-1 atoms 
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FIG. 3: Calculated maximum vibrational frequency (u), Sec¬ 
ond difference of the energy A 2 E{n), absolute binding energy 
per atom Ei,, average coordination number of atoms Wc, and 
the average bond length d of the pure, mono-doped, bi-doped 
Ag and Cu clusters as a function of size. The empty circles 
show the measured binding energy of the pure Cu clusters 

M. 


clusters. This parameter is defined as follows: 

A2E{n) = Etot{n + 1 ) + Etotiji - 1 ) - 2Etot{n) ( 2 ) 

Etot{n) is the minimized energy of the M„_mPdm clus¬ 
ter (M = Ag, Cu) and Etotin + 1) {Etot{n — 1)) is the 
minimized energy of the cluster with one more (less) M 
atom. The calculated second-order energy difference of 
the pure and doped Cu and Ag clusters is plotted in Fig. 
[31 It is seen that pure clusters show strong peaks at 
n = 2, 8 which indicate the highest relative stability of 
these systems and propose n = 2,8 as the first magic 
numbers of small pure Ag and Cu clusters, in agreement 
with the experimental observations and electronic 
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FIG. 4: Calculated HOMO-LUMO gap, Ionization potential 
(IP), and electron affinity (EA) of the pure, mono-doped, and 
bi-doped Ag (left column) and Cu clusters (right column) 
after the GW correction. The empty circles show the available 
measured data for pure clusters, collected from references [H- 

m. 


shell jellium model [13, • The predicted magic num¬ 

bers for mono-doped clusters are 3 and 7, while in the 
case of bi-doped clusters, the main peaks are observed at 
n = 3,4, 8. 

Furthermore, a clear odd-even oscillation is seen in the 
pure and mono-doped clusters, although bi-doped clus¬ 
ters do not display such clear oscillations. It is worth¬ 
while to mention that electron pairing in HOMO can ex¬ 
plain this oscillations and the more stability of the even 
pure and odd mono-doped clusters. The valence elec¬ 
trons of the pure M„ and mono-doped M„_iPd clusters 
are mainly composed of the M s electrons, hence clusters 
with odd number of M atoms have one unpaired elec¬ 
tron and consequently exhibit lower stability. The pure 
and mono-doped clusters with even number of M atoms 
have a close electronic shell and hence are more stable. 
On the other hand, the behaviour of the A 2 E curve of 
the bi-doped clusters is more complicated and needs to 
further arguments. We attribute this more complicated 
behaviour to the effective contribution of d electrons to 
the valence shell of the bi-doped clusters. As it was ar¬ 
gued in section III, the direct bonding between Pd atoms 
in the bi-doped clusters give rises to effective promotion 
of the Pd d electrons to the valence shell of the system. 
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FIG. 5: Calculated HOMO-LUMO gap of the pure Ag and 
Gu clusters before and after the GW correction along with the 
available experimental data, collected from references 
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VI. ELECTRONIC PROPERTIES 


Since the single particle Kohn-Sham eigenvalues have 
no clear physical meaning, the reliable energy level of 
the highest occupied (HOMO) and the lowest unoccu¬ 
pied molecular orbital (LUMO) of the clusters are deter¬ 
mined after applying the many-body GW perturbation 
theory [ 2 ^. The GW corrected HOMO and LUMO en¬ 
ergies are expected to be comparable to the experimen¬ 
tal ionization potential (IP) and electron affinity (EA) 
of the clusters. These values along with the computed 
HOMO-LUMO gaps are presented in Fig. 0] The avail¬ 
able measured data for pure clusters are also presented 
in this figure. The calculated values of IP and Gap ex¬ 
hibit good agreement with the experimental data which 
is mainly due to the GW correction, while the calculated 
EA values in some cases are considerably overestimated. 
In order to show the effect of this many body correc¬ 
tion, we have compared in Fig. [S]the experimental values 
of the HOMO-LUMO gap of the pure clusters with the 
calculated ones before and after application of the GW 
correction. It is seen that the Kohn-Sham gaps are signif¬ 
icantly underestimated while the GW corrected gaps are 
generally much closer to experiment. The observed odd- 
even oscillations in the IP and EA of the pure clusters 
are attributed to the electron pairing in these systems, 
discussed in the previous section. The lower value of EA 
for even clusters indicates less tendency of these closed 
shell systems to receive an extra electron, compared with 
the open shell clusters. A consistent reversed trend is 
visible in the obtained ionization potentials. 

In the same way, the ionization potential and electron 
affinity of the most stable structure of the doped clusters 
were calculated and displayed in Eig.Sl The higher value 
of IP for the odd M„_iPd clusters indicates their more 
chemical stability compared to the neighboring even clus¬ 
ters. In the bi-doped clusters, the IP and EA curves do 
not show clear oscillations with respect to cluster sizes, 
which is likely due to the effective contribution of d elec¬ 
trons to the Pd-Pd bonding in these clusters. The trend 
of the calculated HOMO-LUMO gap is generally con¬ 
sistent with the A 2 E, indicating that the systems with 
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higher gap are generally more stable. It is consistent 
with the chemical intuition that a higher value of the 
HOMO-LUMO gap indicates less chemical activity and 
hence more stability. Considering the spin polarized elec¬ 
tronic structure of the most stable isomers, we found out 
that consistent with the simple spin pairing rule, all stud¬ 
ied clusters prefer the lowest possible spin multiplicity. 
Hence, the pure and doped clusters with even number of 
Ag and Cu atoms are nonmagnetic while others exhibit 
a total spin moment of 1 /i^. 

VII. CONCLUSIONS 

In this paper, spin-polarized all-electron calculations 
were performed to investigate structural and electronic 
properties of the pure and Pd doped Cu„ and Ag„ clus¬ 
ters (n < 9). It was argued that a 2D to 3D structural 
cross over occurs in the size of seven of the pure clusters. 
The magic numbers of the pure clusters were found to be 
2 and 8. The obtained results indicate that the doped Pd 
atoms generally prefer the high coordinated sites of the 
pure clusters. We found that the magic numbers of the 
mono-doped clusters are 3 and 7 while bi-doped clusters 


exhibit magic numbers 3, 4, and 8. It was argued that 
doping with one Pd atom slightly reduces the onset of 
the 2D to 3D structural crossover of the pure Ag clus¬ 
ters, while doping with two Pd atoms strongly enhances 
stability of the 3D isomers and completely destroys this 
cross over. Moreover, it was observed that the bond stiff¬ 
ness of the Cu clusters are decreased after doping with 
Pd atoms while their average bond strength have lower 
sensitivity to doping. On the other hand, Pd doping in¬ 
creases the bond stiffness and bond strength of the Ag 
clusters. It was discussed that electron pairing in the 
HOMO, induces a significant odd-even fluctuation in the 
ionization potential, electron affinity, and HOMO-LUMO 
gap of the pure Ag and Cu clusters, while doping with 
Pd atoms considerably weakens this fluctuation. 
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